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We study the GOY shell model simulating the cascade processes of turbulent flow. The model has 
two inviscid invariants governing the dynamical behavior. Depending on the choice of interaction 
coefficients, or coupling parameters, the two invariants are either both positive definite, analogous 
to energy and enstrophy of 2D flow, or only one is positive definite and the other not, analogous to 
energy and helicity of 3D flow. In the 2D like model the dynamics depend on the spectral ratio of 
enstrophy to energy. That ratio depends on wave- number as k a . The enstrophy transfer through 
the inertial sub-range can be described as a forward cascade for a < 2 and a diffusion in a statistical 
equilibrium for a > 2. The a = 2 case, corresponding to 2D turbulence, is a borderline between the 
two descriptions. The difference can be understood in terms of the ratio of typical timescales in the 
inertial sub-range and in the viscous sub-range. The multi fractality of the enstrophy dissipation 
also depends on the parameter a, and seems to be related to the ratio of typical timescales of the 
' different shell velocities. 

<n : 

I I. INTRODUCTION 

> ; 

The standard Kolmogorov A; -5 / 3 scaling law for energy cascading in 3D turbulence and the corresponding k~ 3 
scaling law for enstrophy cascade in 2D turbulence is still debated. Direct numerical calculations of the full Navier- 
Stokes equation is by and large still impossible for high Reynolds number (> 100 — 200) flows. However, the cascading 
t-H ' mechanisms and its multi fractal nature can be analyzed in reduced wave-number models for very high Reynold 
numbers with high accuracy. In this paper we investigate the GOY shell model |l],|| where the spectral velocity or 
vorticity is represented by one complex variable for each shell evenly spaced in log(k) in spectral space. For this type 
of model the Kolmogorov scaling arguments can be applied as for real flow regardless of how realistically they mimic 
the dynamics of the Navier-Stokes equation. The scaling behavior of the fields depends on the inviscid invariants of 
the model. In the simple model we are able to control which symmetries and conserved integrals of the dynamics that 
are present in the inviscid and force- free limit. In the models we interpret as simulating 3D turbulence there are 2 
inviscid invariants, similar to energy and helicity |3j, of which the first is positive definite and the second is not. For 
the models we interpret as 2D turbulence the 2 inviscid invariants, similar to energy and enstrophy, are both positive 
definite. We will mainly be concerned with an investigation of the 2D like models. The specific parameter choice 
previously assigned to simulating 2D turbulence are such that the GOY model does not show enstrophy cascading 
but rather a statistical equilibrium where the enstrophy is transported through the inertial sub-range by diffusion 
j|. We show that this is a borderline case for which, on one side, the model behaves as a cascade model and, on 
the other side, it behaves as a statistical equilibrium model, where the enstrophy spectrum is characterized by a 
simple equipartitioning among the degrees of freedom of the model. The difference in behavior is connected with 
the different typical timescales of the shell velocities as function of shell number. This probabily also influences the 
(non-universal) multi fractal behavior of the shell velocities. If timescales in the viscous sub-range are not smaller 
than in the beginning of the inertial sub-range, the low wave-number end, the model does not have a multi fractal 
spectrum. 
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II. THE GOY MODEL 



The GOY model is a simplified analogy to the spectral Navier-Stokes equation for turbulence. The spectral domain 
is represented as shells, each of which is defined by a wavenumber k n = k$\ n , where A is a scaling parameter defining 
the shell spacing; in our calculations we use the standard value A = 2. The reduced phase space enables us to cover a 
large range of wavenumbers, corresponding to large Reynolds numbers. We have 2N degrees of freedom, where TV is 
the number of shells, namely the generalized complex shell velocities or vorticities, u n for n = l,N. The dynamical 
equation for the shell velocities is, 
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ii n — ik n (au n+2 u n+1 + -^u n+1 u n _ 1 + ~^u n _- i u n _ 2 ) ~~ vk^Un — v k n P2 u n + f n , (1) 

where the first term represents the non-linear wave interaction or advection, the second term is the dissipation, the 
third term is a drag term, specific to the 2D case, and the fourth term is the forcing. Throughout this paper we 
use pi = p 2 = 2. We will for convenience set a = ko = 1, which can be done in (Q) by a rescaling of time and the 
units in k space. A real form of the GOY model, as originally proposed by Gledzer [jj], can be obtained trivially 
by having purely imaginary velocities and forcing. The GOY model in its original real form contains no information 
about phases between waves, thus there cannot be assigned a flow field in real space to the spectral field. The complex 
form of the GOY model and extensions in which there are more shell variables in each shell introduce some degrees 
of freedom, which could be thought of as representing the phases among waves. However, it seems as if these models 
do not behave differently from the real form of the model in regard to the conclusions in the following . The key 
issue for the behavior of the model is the symmetries and conservation laws obeyed by the model. 



A. Conservation laws 



The GOY model has two conserved integrals, in the case of no forcing and no dissipation (y = f = 0). We denote 
the two conserved integrals by 

N N N 

= J2 E n 2 = lJ2 C< 2 Kl 2 = ~ £ x nai ' 2 K| 2 (2) 

n—1 n— 1 n— 1 

By setting E 1,2 =0 and using ii n from (|l|) we get 

f 1 + bzi + cz\ = 
\ 1 + bz 2 + czl = 0, 



(3) 

where the roots zi^ = A" 12 are the generators of the conserved integrals. In the case of negative values of z we can 
use the complex formulation, a = (log\z \ + inj/logX. The parameters (b, c) are determined from (^) as 

b = -{zi + z 2 )/z 1 z 2 
c = l/zxz 2 . 



(4) 

In the (6, c) parameter plane the curve c = b 2 /4 represents models with only one conserved integral, see figure 1. 
Above the parabula the generators are complex conjugates, and below they are real and different. Any conserved 
integral represented by a real nonzero generator z defines a line in the (b, c) parameter plane, which is tangent to the 
parabula in the point (b, c) — (— 2/z, 1/z 2 ). The rest of our analysis we will focus on the line defined by z\ = 1. The 
conserved integral, 

1 - 

£ 1 = 2Ei u «i 2 ' ( 5 ) 

n=l 

is the usual definition of the energy for the GOY model ||. The parameters are then determined byl + & + c = 0, 
which with the definitions b — — e and c — — (1 — e) agree with the notation of ref. ||. The generator of the other 
conserved integral is from (Q) given as, 

* - ^j. (0) 
For e < 1 the second conserved integral is not positive definite and is of the form, 
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E 2 =H=±^2(-ir\z 2 \ n \u n \ 2 , (7) 

n=l 

which can be interpreted as a generalized helicity. For e = 1/2, Z2 = — 2 = — A the model is the usual 3D shell model 
and H is the helicity as defined in ref. By choosing A such that A = 1/(1 — e) we get E 2 = ^(—l) n X n \u n \ 2 . This 
form was argued in ref. || to be the proper form for the helicity. In this paper we will alternatively use the definition 
(fjj) for the helicity. 

For e > 1 the second conserved integral is positive definite and of the form, 



E 2 = Z = 



1 N 

^E^KI 2 , (8) 



=i 



which can be interpreted as a generalized enstrophy. For e = 5/4, Z2 = 4 = A 2 the model is the usual 2D shell 
model and Z is the enstrophy as defined in ref. Q . The sign of c, which is the interaction coefficient for the smaller 
wavenumbers, changes when going from the 3D - to the 2D case. This could be related to the different role of backward 
cascading in the two cases. To see this, consider the non-linear transfer of E l through the triade interaction between 
shells, n — 1, n, n+ 1. This is simply given by, 



E n -1 


- h ai A 


K 


- hy h ai A 


E n+1 


_ 2lo, a 



(9) 
with 

A„ = k n -xlm(u n -iu n u n+ i). (10) 

The detailed conservation of E l in the triade interaction is reflected in the identity, 1 + bzi + cz 2 = 0. Using (^) 
and (p) we have for the exchange of energy, E 1 , with a\ = 0; 



E n-1 - &n 

E n = -eA„ 

E n+1 = (e - 1)A„ (11) 



and for helicity/enstrophy, E 2 , with a.2 = a; 

E? — h a A 

S 2 = -(e/(e-l))fc^ 1 A„ 
^ + i = (l/(e-l))fc°_ 1 A„. (12) 

We have e < 1 for 3D like models and e > 1 for 2D like models, the two situations are depicted in figure 2, where 
the thickness of the arrows symbolize the relative sizes of the exchanges in the cases of e = 1/2 and e = 5/4. 



B. Scaling and inertial range. 

The inertial sub-range is defined as the range of shells where the forcing and the dissipation are negligible in 
comparison with the non-linear interactions among shells. Since we apply the forcing at the small shell numbers and 
the dissipation at the large shell numbers, the inertial range (of forward cascade) is characterized by the constant 
cascade of one of the conserved quantities. The classical Kolmogorov scaling analysis can then be applied to the 
inertial range. There is, however, in the shell model, long range influences of the dissipation and forcing into the 
inertial subrange. This is an artifact of the modulus 3 symmetry, see equation (pp|), and the truncation of the shell 
model which is not expected to represent any reality. These features are treated in great detail in ref. |^]. Denoting 
7712 as the average dissipation of E 1,2 , this is then also the amount of E 1 ' 2 cascaded through the inertial range. The 
spectrum of E ' does then, by the Kolmogorov hypothesis, only depend on k and 7/1,2- From dimensional analysis 
we have, [ku] = s _1 , [771,2] = [E 1 ' 2 ^' 1 , [E 1 - 2 ] = [k a ^u 2 ] = [k] a - 2 S - 2 , and we get, 
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E 1 ' 2 ~ ^fe^.- 2 )/ 3 . (13) 

For the generalized velocity, u, we then get the " Kolmogorov- scaling", 

|u| ~ ry^V^^H 1 )/ 3 . (14) 

The non-linear cascade, or flux, of the conserved quantities defined by zi, 2 through shell number n can be expressed 
directly as, 



n 1 n 

II*' 2 = E h2 (m) = ~£ <2««™ + c.c.) 



2 

m— 1 



z™ 2 (-A„/z 24 + A„+i). (15) 



In the inertial range the cascade is constant, II*' 2 = n n ', so from (15) we get following ref. 



+i 

ziz 2 A n+2 - [zi + z 2 )A„+i + A„ = => 

q n +ziz 2 /q n +i = zi +%2 (16) 

where we have defined 

q n = A n /A n+1 . (17) 
The inertial range scaling requires q n = q n +i = q to be independent of n. Solving (|l^) for q and using (|l7|) gives, 

{Zi => u„ ~ fen ( - Q1+1 ^ 3 Kolmogorov for i? 1 
^2 => % ~ fe„ a2+1 ^' 3 Kolmogorov for i? 2 . 

Inserting this into ( |T^ ) gives for the cascade of -E 1 in the two solutions, 



! J 1 — z 2 /zi Kolmogorov for i? 1 
" < fluxless for B 1 , 



(18) 



and correspondingly for E 2 , 



n 2 



fluxless for E 2 . 

1 — zi/z 2 Kolmogorov for E 2 



These are the two scaling fixed points for the model. The Kolmogorov fixed point for the first conserved integral 
corresponds to the fluxless fixed point for the other conserved integral and visa versa. This is of course reflected in 
the fact that ( [L6| ) is symmetric in the indices 1 and 2. That these points in phase space are fixed points, in the case 
of no forcing and dissipation, is trivial, since n„ = n n+ i => E n+ \ = => u n +i = 0. It should be noted that the 
Kolmogorov fixed point, 

u-k-^+V/ 3 , (19) 

obtained from this analysis is in agreement with the dimensional analysis (fl3|). 

The scaling fixed points can be obtained directly from the dynamical equation as well. For u n ~ k~ 7 g(n) = 
\~ ni g(n), where g(n + 3) = g(n) is any period 3 function, we get by inserting into (0) with a = 1, 

g(n - l)g(n)g(n + ljA^-^+^tl + &A 37 " 1 + c(A 37 ~ 1 ) 2 ) = (20) 

and the generators reemerge, zi. 2 = A" 12 = A 7l,2_ , giving the Kolmogorov fixed points for the two conserved 
integrals, 71,2 = (0:1,2 + l)/3. 

The period 3 symmetry seems to have little implications for the numerical integrations of the model, except perhaps 
in accurately determining the structure function. 

The stability of the fixed point for energy cascade in the 3D case, e = 1/2, is characterized by few unstable directions, 
where the corresponding eigenmodes mainly projects onto the high shell numbers, and a large number of marginally 
stable directions which mainly projects onto the inertial range. This also holds in the case with forcing and dissipation 
^ . In the case where dissipation and forcing are applied for some values of e the Kolmogorov fixed point can become 
stable. Biferale et al. ra] show that there is a transition in the GOY model, for v = 10 -6 and / = 5 x 10~ 3 x (1 +%), 
as a function of e from stable fixed point (e < 0.38..), through Hopf bifurcations and via a Ruelle-Takens scenario to 
chaotic dynamics (e > 0.39..). 
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III. FORWARD AND BACKWARD CASCADES 



Until this point we have not specified which of the two conserved quantities will cascade. Assume, in the chaotic 
regime where the Kolmogorov fixed points are unstable, that there is, on average, an input of the same size of the 
two quantities, E 1 and E 2 , at the forcing scale, this can of course always be done by a simple rescaling of one of 
the quantities. If is a shell number at the beginning of the viscous subrange, we have that u^ d /k]^ d w ^,and the 
dissipation, D\ of the conserved quantity, E z , can be estimated as 

Z^~^ +2 W| 2 . (21) 

The ratio of dissipation of E 1 and E 2 scales with k^ d as D 1 /D 2 <~ fc^ - " 2 , so that, in the limit Re — > oo when 
a.\ < ct2, there will be no dissipation in the viscous sub-range of E 1 where E 2 is dissipated. Therefore, a forward 
cascade of E 1 is prohibited and we should expect a forward cascade of E 2 . For the backward cascade the situation is 
reversed, so we should expect a backward cascade of E 1 . 

The situation is completely different in the 2D like and the 3D like cases. In the 3D like models E 2 is not positive 
definite, E 2 (helicity) is generated also in the viscous sub-range and for the usual GOY model we do not see a forward 
cascade of helicity, see, however, ref. 0. This is in agreement with the observed A: -5 / 3 energy spectrum observed in 
real 3D turbulence corresponding to the forward cascade of energy. In the 2D case we observe the direct cascade of 
enstrophy, while the inverse cascade of energy is still debated. In the rest of this paper we will concentrate on 2D like 
models where we will implicitly think of E 1 = E, with ot\ — 0, as the energy and E 2 — Z , with a.2 — a > 0, as the 
enstrophy. With regard to the inverse cascade of energy one must bare in mind that in 2D turbulence the dynamics 
involved is probably related to the generation of large scale coherent structures, vortices, and vortex interactions. 
Vortices are localized spatially, thus dclocalized in spectral space. This is in agreement with the estimate that 2D is 
marginally delocalized in spectral space . In the GOY model there is no spatial structure and the interactions are 
local in spectral space. The model is therefore probably not capable of showing a realistic inverse energy cascade. We 
will thus only consider the forward cascade in this paper. Figure 3 shows the scaling in the inertial sub-range of the 
model with e = 5/4 corresponding to a = 2. The cascades of the enstrophy and energy are shown in figure 4. It is 
seen that enstrophy is forward cascaded while energy is not. 



IV. STATISTICAL DESCRIPTION OF THE MODEL 



In a statistical equilibrium of an ergodic dynamical system we will have a probability distribution among the (finite) 
degrees of freedom, assuming an ultraviolet cutoff, of the form, P t ~ exp(—BEl — AE 2 ), where E 1 and E 2 are the 
conserved quantified, energy and enstrophy. Thus, the temporal mean of any quantity, which is a function of the shell 
velocities is given as 

9 = J Y\_ diiig(u\, ...,un) exp(-BEl - AE 2 )/ J \{ du t exp(-BEl - AE 2 ). (22) 

i i 

A and B are Lagrange multipliers, reflecting the conservation of energy and enstrophy when maximizing the entropy 
of the system, corresponding to inverse temperatures, denoted as inverse "energy-" and " enstrophy-temperatures" 
. The shell velocities themselves will in this description be independent and gaussian distributed variables with 
standard deviation tr(uj) = l/(2(_Bfc" 1 + Ak" 2 )). The average values of the energy and enstrophy becomes, 

e} = k?W =( B + Ak^-^y 1 

E 2 = k^\^=(Bk^- a2 +A)-\ (23) 

For k — y we will have equipartitioning of energy, k" 1 \ui\ 2 = B^ 1 and the scaling ~ k i ai ^ 2 and for the other 
branch, k — > oo, we will have equipartitioning of enstrophy /c" 2 |wi| 2 = A^ 1 and the scaling \m\ ~ k i " 2 / 2 . I n the 
case of no forcing and no viscosity the equilibrium will depend on the ratio A/B between the initial temperatures 
A^ 1 , B^ 1 . To illustrate this we ran the model without forcing and viscosity but with 2 different initial spectral slopes 
of the velocity fields, the larger the slope the higher the ratio of the energy temperature to the enstrophy temperature. 
Figure 5 shows the equilibrium spectra for e = 5/4, v = f = 0, in the cases of initial slopes -1, -0.8. The full lines are 
the equilibrium distribution given by (^3| ) for A/B — 10 2 and A/B — 10~ 2 respectively. 
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V. DISTINGUISHING CASCADE FROM STATISTICAL EQUILIBRIUM 



For the forward enstrophy cascade the spectral slope is — (a + l)/3 and the enstrophy equipartitioning branch has 
spectral slope —a/2. Thus for the 2D case where a = 2 we cannot distinguish between statistical (quasi-) equilibrium 
and cascading. This was pointed out by Aurell et al. Q and it was argued that the model can be described as being 
in statistical quasi-equilibrium with the enstrophy transfer described as a simple diffusion rather than an enstrophy 
cascade. This coinciding scaling is a caviate of the GOY model not present in the real 2D flow where the statistical 
equilibrium energy spectrum scales as fc _1 and the cascade energy spectrum scales as fc -3 . For other values of a the 
scaling of the two cases arc different, see figure 6. This figure represents the main message of this paper. First axis is 
the parameter e, along the line shown in fig. f , defining the spectral ratio between the two inviscid invariants. Second 
axis is the scaling exponent 7. The horizontal dashed line 7 = 1/3 is the Kolmogorov scaling exponent for energy 
cascade. The full curve is the scaling exponent for the enstrophy cascade, and the dotted curve corresponds to the 
enstrophy equipartitioning. 

All the 3D like models (asterisks in figure 6) are near energy cascade scaling (dashed line). Statistical equilibrium 
corresponds to the line 7 = 0. The bold line piece, < e < 0.39..., represents parameter values where the Kolmogorov 
fixed point is stable [||. The scaling for e > 0.39... is slightly steeper than the Kolmogorov scaling, which is attributed 
to intermittency corrections originating from the viscous dissipation |l2] | . ft seems as if there is a slight trend showing 
increasing spectral slopes for increasing e. 

For the 2D like models the scaling slope is also everywhere on or slightly above both the cascade - and the equilibrium 
slopes (diamonds in the figure). The classical argument for a cascade is that given an initial state with enstrophy 
concentrated at the low wave-number end of the spectrum, the enstrophy will flow into the high wave-numbers in order 
to establish statistical equilibrium. The ultra-violet catastrophe is then prevented by the dissipation in the viscous 
sub-range. Therefore, we cannot have a non-equilibrium distribution with more enstrophy in the high wave-number 
part of the spectrum than prescribed by statistical equilibrium since enstrophy in that case would flow from high - to 
low wave-numbers. This means that the spectral slope in the inertial sub-range always is above the slope corresponding 
to equilibrium (dotted line in figure 6). Consequently, the 2D model with e = 5/4 separates two regimes, 1 < e < 5/4 
where enstrophy equilibrium is achieved and 5/4 < e < 2 where the enstrophy is cascaded through the inertial range. 

In figure 7 the spectra and the cascades are shown for different values of e. The model was run with 50 shells and 
forcing on shell number 15 for 2 x 10 4 time units and averaged. Even then there are large fluctuations in the cascades 
not reflected in the spectra. The large differences in the absolut values for the cascades, Pi, is a reflection of the 
scaling relation (|2l|). 

We interpret the peaks around the forcing scale for e — 11/10 as statistical fluctuation and the model shows no 
cascade. For e > 5/4 we see an enstrophy cascade and what seems to be an inverse energy cascade. However, we 
must stress that we do not see a second scaling regime for small n corresponding the inverse cascade. Note that for 
e = 2 energy and enstrophy are identical and we have only one inviscid invariant. So if a regime of inverse energy 
cascading existed in parameter space near e — 2 the scaling exponents will be almost identical and coincide at e — 2. 

The two regimes corresponding to equipartitioning and cascade can be understood in terms of timescales for 
the dynamics of the shell velocities. A rough estimate of the timescales for a given shell n, is from (|) given as 
T n ~ (k n u n ) ~ K^ 1 . Again e = 5/4, corresponding to 7 = 1, becomes marginal where the timescale is independent 
of shell number. For e < 5/4 the timescale grows with n and the fast timescales for small n can equilibrate enstrophy 
among the degrees of freedom of the system before the dissipation, at the "slow" shells, has time to be active. 
Therefore these models exhibits statistical equilibrium. For e > 5/4 the situation is reversed and the models exhibits 
enstrophy cascades. Time evolutions of the shell velocities are shown in figure 8, where the left columns show the 
evolution of a shell in the beginning of the inertial subrange and the right columns show the evolution of a shell at 
the end of the inertial subrange. This timescale scaling might also explain why no inverse cascade branch has been 
seen in the GOY model. The timescales at the small wave-number end of the spectrum, with the dissipation or drag 
range for inverse cascade, is long in comparison with the timescales of the inertial range of inverse cascade. Therefore 
a statistical equilibrium will have time to form. The analysis suggests that parameter choices e > 5/4 might be more 
realistic than e = 5/4 for mimicing enstrophy cascade in real 2D turbulence. 

VI. INTERMITTENCY CORRECTIONS 

The numerical result that the inertial range scaling has a slope slightly higher than the K 41 prediction, is not 
fully understood. This is attributed to intermittency corrections originating from the dissipation of enstrophy in the 
viscous subrange. 
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The evolution of the shell velocities in the viscous sub-range is intermittent for e > 5/4, where the PDF's are 
non-gaussian, while the PDF's for e = 5/4 are gaussian in both ends of the inertial sub-range, see figure 9. The 
deviation from the Kolmogorov scaling is expressed through the structure function, [jl2). The structure function 
is defined through the scaling of the moments of the shell velocities; 

~ (24) 

where 5£(q) is the deviation from Kolmogorov scaling. The structure function, £((?), ancl ^C(s) f° r e — 
11/10,5/4,3/2,7/4,2 are shown in figure 10. For e > 5/4 there are intermittency corrections to the scaling in 
agreement with what the PDF's show. 

We know of no analytic way to predict the intermittency corrections from the dynamical equation. Our numerical 
calculations suggest that the intermittency corrections are connected with the differences in typical timescales from 
the beginning of the inertial sub-range, where the model is forced, to the viscous sub-range. The ratio of timescales 
between the dissipation scale and the forcing scale can be estimated by; T v /Tf ~ \ AN ( 1 +i) j where AN is the number 
of shells beween the two. Figure 11 (a) shows the numerical values of #C(10) as a function of e and figure 11 (b) shows 
log2(T v /Tf) as a function of e. The vertical line indicates the crossover between statistical equilibrium and cascading. 
We must stress that caution should be taken upon drawing conclusions from this since the authors have no physical 
explanation of the appearent relationship. 



VII. SUMMARY 



The GOY shell model has two inviscid invariants, which govern the behavior of the model. In the 2D like case 
these corresponds to the energy and the enstrophy of 2D turbulent flow. In the model we can change the interaction 
coefficient, e, and tune the spectral ratio of enstrophy to energy, Z n /E n = k™. For a > 2 we can describe the 
dynamics as being in statistical equilibrium with two scaling regimes corresponding to equipartitioning of energy and 
enstrophy respectively. The reason for the equipartitioning of enstrophy in the inertial range (of forward cascading of 
enstrophy) is that the typical timescales, corresponding to eddy turnover times, are growing with shell number, thus 
the timescale of viscous dissipation is large in comparison with the timescales of non-linear transfer. Thus, this choice 
of interaction coefficient is completely unrealistic for mimicing cascades in 2D turbulence. For a < 2 the model shows 
forward cascading of enstrophy, but we have not identified a backward cascade of energy. The usual choice e = 5/4, 
a = 2 is a borderline and we suggest that a < 2 in respect to mimicing enstrophy cascade might be more realistic. We 
observe that the dynamics becomes more intermittent when a < 2, in the sense that the structure function deviates 
more and more from the Kolmogorov prediction. For e = 2 we have a = 0, thus energy and enstrophy degenerates 
into only one inviscid invariant, this point could then be interpreted as a model of 3D turbulence. However, as is 
seen from (^), in this case the fluxless fixed point is the one surviving, but as is seen in figure 7, bottom panels, 
this model also shows cascading. This choice for 3D turbulence model could shed some light on the dispute of the 
second inviscid invariant (helicity) being important || or not [fl3| for the deviations from Kolmogorov theory, work is 
in progress on this point. 
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FIG. 1. The (b, c) parameter space. The line indicates where the one conserved integral is the energy, defined as 

E=l/2j2Wn\ 2 
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2D: 
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FIG. 2. The shell triad interactions. Arrows indicates transfer of energy, E 1 , and enstrophy, E 2 , for the 2D case, e = 5/4, 
and energy, E 1 , and helicity, E 2 , for 3D, e = 1/2, case. The thickness of the arrows indicates the strength of the transfer. 
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FIG. 3. The mean value of the shell velocities as a function of shell number on a logarithmic scale (base A), for the 2D case, 
e = 5/4, fc = A~ 4 , A = 2, n = 30, v = 1CT 16 , /„ = 5 x 1CT 3 x (1 + i)<5„, 4 . The model was run for 4.2 x 10 4 time-units. 
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FIG. 4. The mean values of, (a), the energy flux, III, and (b) the enstrophy flux, II2. 



12 




FIG. 5. Same as figure 3, but with n = 30, v = / = 0. Diamonds corresponds to an initial spectral slope of -1.0, that is a 
high value of A/B. The corresponding curve is statistical equilibrium distribution for A/B = 10 2 . Triangles corresponds to an 
initial spectral slope of -0.8, that is a lower value of A/B. The curve is statistical equilibrium distribution for A/B = 10~ 2 . 
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FIG. 6. The spectral slope 7 as function of e. The horizontal dashed line is the Kolmogorov scaling for 
energy cascade. The full curve is the scaling exponent for the enstrophy (or in 3D like case helicity) cascade, 
and the dotted curve corresponds to the enstrophy (helicity) equipartitioning. All the 3D like models shows en- 
ergy cascade (equipartitioning corresponds to the line 7 = 0). The bold line piece, < e < 0.39..., rep- 
resents parameter values where the Kolmogorov fixed point is stable |J. The diamonds are model run with 
n = 50,fc = A -4 , A = 2,/„ = 5xl(T 4 x(l+i)<5 n ,i 5 , (e = 11/10, 1^ = 5xl(T 27 y = 100), (e = 117/100, v = 5xl0~ 27 ,^' = 100), (e = 5/4, v = I 
The stars are model runs with n = 19, k = A" 4 , A = 2, f n = 10~ 4 x (1 + i)5„ A , v = 10~ B , v' = 0, e = 1/2, 6/10, 7/10. The 3D 
like models shows Kolmogorov scaling, with deviations due to intermittency corrections, and energy cascading. The 2D like 
models shows a cross-over at e = 5/4 between statistical equilibrium, 1 < e < 5/4, and enstrophy cascading, 5/4 < e < 2. 
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FIG. 8. Time evolution of shell velocities in the beginning and the end of the inertial subrange. The typical timescale of 
shell n scales as T n ~ (fc n |w n |) _1 ~ k^ 1 . Note that for e = 5/4 the timescale is the same for all shells. 
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FIG. 10. The deviation of the structure function from Kolmogorov scaling for e = 11/10, 5/4, 3/2, 7/4, 2. 
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FIG. 11. (a) Numerical values of <5C(10) as a function of e, the error bars represents one standard deviation, (b) Ratio of 
typical timescales between dissipation scale and forcing scale, log^iTvjTf) as a function of e. The vertical line indicates the 
crossover between statistical equilibrium and cascading. 
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